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1  Introduction 


The  collapse  of  vapor  cavities  in  the  vicinity  of  rigid  walls  has  been  the  subject  of  a  number 
of  theoretical,  numerical  and  experimental  investigations  [1,2, 3, 4, 5, 6].  The  collapse  is  char¬ 
acterized  by  a  jet  which  is  formed  at  the  side  of  the  cavity  farthest  from  the  rigid  wall  and 
is  directed  toward  the  wall.  It  is  believed  that  this  jet  enhances  the  erosion  of  solid  material 
from  the  wall.  Several  investigators  have  reported  that  laboratory  cavitation  erosion  tests 
of  relatively  rigid  surfaces  coated  with  soft  elastomeric  materials  show  increased  resistance 
to  erosion  [7,8,9].  To  explore  the  mechanism  for  this  phenomenon  Gibson  and  Blake  [10] 
photographed  the  collapse  of  a  cavity  in  the  vicinity  of  walls  covered  by  elastomeric  coatings 
of  several  compositions.  They  found  that  under  some  circumstances  the  the  collapse  wcis 
modified  such  that  a  jet  was  formed  on  the  side  of  the  cavity  closest  to  the  wall  and  was 
directed  away  from  the  wall.  It  is  theorized  that  this  redirection  of  the  jet  protects  the  wall 
from  erosion.  To  analyze  this  effect  theoretically  they  devised  a  simple  model.  First,  using 
data  from  the  collapse  of  a  spherical  cavity  in  an  infinite  fluid,  they  computed  the  average 
pressure  force  and  average  displacement  on  an  imaginary  surface  at  the  position  where  the 
wall  would  be  located  in  the  interactive  case.  Then  using  a  one-dimensional  spring-mass- 
damper  model  of  the  coating,  they  found  coating  designs  whose  impedance  matched  the 
average  force  and  average  displacement  characteristics  in  the  spherical  collapse.  In  this  way 
they  attempted  to  find  the  coating  properties  that  would  make  the  collapse  spherical  when 
the  cavity  wais  adjacent  to  the  compliant  wall.  It  was  theorized  that  further  changes  in  the 
wall  properties  would  direct  the  jet  away  from  the  wall  as  occurs  when  the  cavity  collapses 
near  a  free  surface.  This  analysis  was  a  good  first  attempt  and  indicated  a  tentative  benefi¬ 
cial  parameter  range  for  the  coatings. 

The  computational  models  for  the  collapse  of  a  cavity  in  the  vicinity  of  a  rigid  wall  [1,2,5] 
or  a  free  surface  [11,12]  have  become  quite  advanced.  In  these  models  the  fluid  motion  is  as¬ 
sumed  to  be  inviscid  and  incompressible  and  the  fluid  velocity  is  obtained  from  the  gradient 
of  a  scalar  Eulerian  velocity  potential  which  satisfies  Laplace’s  equation.  Fluid  particles  on 
the  free  surfaces  are  tracked  in  time  and  Bernoulli’s  equation  is  used  to  obtain  the  velocity 
potential  at  the  position  of  these  particles.  At  each  time  step  an  integral  equation  is  solved 
to  obtain  the  value  of  the  fluid  velocity  in  the  direction  normal  to  the  surfaces.  The  rigid 
surface  is  usually  simulated  with  an  image  cavity. 


In  the  present  paper,  in  order  to  compute  the  collapse  of  a  cavity  near  a  compliant  wall, 
the  above  model  for  the  fluid  motion  has  been  coupled  to  a  simple  model  of  a  compliant 
coating.  Specifically,  the  coating  is  modeled  as  a  membrane  with  a  spring  foundation.  This 
coating  is  characterized  by  its  mass  per  unit  area  (m),  membrane  tension  (T)  and  spring 
stiffness  per  unit  area  (K).  The  coating  is  coupled  to  the  flow  model  through  the  normal 
velocity  and  the  pressure  at  the  flow-coating  interface.  Linearized  boundary  conditions  are 
applied  at  the  location  of  the  undisturbed  membrane  surface.  This  model  is  fully  interactive. 
It  is  show'n  later  in  this  paper  that  there  are  four  dimensionless  parameters  controlling  the 
characteristics  to  the  collapse.  These  are  Rq/ Zo,TmlTo,  B^K fT  and  KRc/AP  where  Rq 
is  the  initial  cavity  radius,  Zq  is  the  initial  distance  of  the  center  of  the  cavity  from  the 
wall,  To  =  Rq^pUPo:,  -  Pq)  is  the  time  scale  for  collapse  of  the  cavity  in  an  infinite  fluid, 
Tm  =  yJm/K  is  the  characteristic  spring-mass  time  scale  for  the  coating,  p  is  the  density 
of  the  fluid,  Pcx,  is  the  pressure  of  the  fluid  far  from  the  cavity,  and  Pq  is  the  pressure  in 
the  cavity  which  is  assumed  constant.  Results  of  the  first  computations  from  the  model  are 
presented  herein.  In  these  calculations,  Rof Zo,TmlTo  and  RlK jT  were  held  constant  at 
1.5, 1.0  and  12.5  respectively  while  KRq/AP  was  varied  from  oo  (a  rigid  wall)  to  12.5.  More 
extensive  calculations  are  underway  at  the  present  time. 

The  remainder  of  the  paper  is  organized  as  follows.  In  Section  2,  a  description  of  the 
mathematical  formulation  of  the  full  interaction  problem  is  given.  This  is  followed  in  Section 
3  by  a  description  of  the  numerical  implementation  of  this  theory.  In  Section  4,  the  results 
are  presented  and  finally,  the  conclusions  of  the  work  are  given  in  Section  5. 

2  Mathematical  Formulation 

A  schematic  showing  the  initial  position  of  the  cavity,  the  compliant  wall  and  the  coordinate 
system  used  in  the  calculations  is  given  in  Figure  1.  A  cylindrical  coordinate  system  is 
used  with  the  z-axis  piercing  the  center  of  the  cavity  and  directed  normal  to  the  undisturbed 
surface  of  the  compliant  wall  which  is  located  at  z  =  0.  In  this  coordinate  system  the  problem 
is  axisymmetric  about  the  z-axis.  The  cavity  is  initially  {t  =  0)  spherical  with  radius  Ro, 
its  maximum  value,  and  its  center  is  located  at  Zq.  For  times  before  the  initial  instant  the 


pressure  everywhere  in  the  fluid  i  Poo  and  the  cavity  surface  is  rigid.  The  internal  pressure 
in  the  cavity  is  assumed  constant,  Pq,  for  all  time.  For  /  >  0  the  boundary  of  the  cavity  is 
a  free  surface  and  the  pressure  in  the  fluid  far  from  the  cavity  is  constant  with  value  Poo- 
The  theoretical  model  for  the  fluid  is  similar  to  the  one  used  by  a  number  of  investigators 
including,  most  recently,  Blake,  Taib,  and  Doherty  [11].  The  fluid  motion  is  assumed  to  be 
incompressible  and  inviscid  and  therefore  satisfies  Laplace’s  equation: 

=  0  (1) 


where  V  is  the  gradient  operator,  ^  is  the  the  velocity  potential  and  u  =  y  ^  where  u  is  the 
fluid  velocity.  On  the  surface  of  the  cavity  the  pressure  in  the  fluid  must  equal  the  pressure 
in  the  cavity,  Pq.  This  can  be  written  as  Bernoulli’s  equation  in  material  derivative  form; 


D4>  ^^n2.Poo-Po 


(2) 


where  Dj Di  is  the  derivative  with  respect  to  time  following  a  fluid  particle.  The  kinematic 
boundary  condition  on  the  surface  of  the  cavity  states  that  material  points  remain  on  the 
surface  of  the  cavitv: 

Dxo 


Dt 


=  V<l> 


(3) 


where  Xc  is  the  position  vector  to  these  material  points. 


The  wall  is  modeled  as  a  membrane  with  tension  T  and  mass  per  unit  area  m  that  is 
supported  by  a  spring  foundation  of  spring  constant  K .  The  equation  describing  the  motion 
of  the  membrane  is 

1  fi  Pin 

(4) 


-  P +  =  -Pdr,  t) 


dP  "  r  dr  ^  dr' 

where  t]  =  »;(r,  t)  is  the  vertical  displacement  of  the  membrane  surface  and  Pc  is  the  pressure 
on  the  membrane.  The  fluid  and  the  membrane  are  coupled  using  linearized  equations  for  the 


pressure  and  velocity  in  the  two  systems.  These  equations  are  satisfied  at  the  undisturbed 
position  of  the  coating  surface  z  —  0: 


Pc(r,f)  =  P(r,0,0 


(5) 

(6) 


For  t  <  0  the  the  springs  are  compressed  uniformly  due  to  the  uniform  pressure  Poo  applied 
by  the  fluid.  The  position  of  the  membrane  at  this  time  is  taken  as  2  =  0. 
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The  time  marching  procedure  to  solve  this  problem  is  as  follows.  Assume  that  at  time 
t  all  dependent  variables  are  known.  The  boundary  conditions  on  the  surface  of  the  cavity, 
Equations  (2)  and  (3),  are  integrated  numerically  to  get  the  position  of  the  surface  of  the 
cavity  and  the  value  of  4>  on  the  cavity  at  time  t  +  At.  The  membrane  equation  can  be  used 
to  obtain  the  value  of  drjfdt  =  d<l>(dn  at  the  membrane  surface  at  t  -f  At.  In  order  to  move 
on  to  the  next  time  step,  t  +  2At,  the  values  of  must  be  known  on  the  cavity  surface  for 
use  in  Equations  (2)  and  (3).  However,  at  this  point  only  the  value  of  d<j>/ds  is  known  (where 
s  is  a  coordinate  along  the  cavity  surface).  Also,  in  order  to  find  the  value  of  d4>ldn  on  the 
membrane  surface  at  1 4-  2At,  the  pressure  must  be  known  on  the  membrane  surface  at  t  +  At 
for  use  in  Equation  (4).  The  pressure  at  t  +  At  can  be  obtained  from  Bernoulli’s  equation  if 
d<l>jdi  and  V*/*  known.  Since  d4>ldn  on  the  membrane  surface  has  been  computed  from 
integrating  the  membrane  equation  in  time  and  (f>  at  time  t  is  known,  only  the  value  of  (j) 
on  the  membrane  at  time  t  +  At  remains  to  be  computed,  then  d<f>lds  and  d<j>ldl  can  be 
computed  from  finite  difference.  Thus,  to  complete  the  problem,  the  values  of  d(l>/dn  on  the 
cavity  surface  and  4>  on  the  membrane  surface  are  obtained  by  solving  Laplace’s  equation  in 
the  form  of  the  integral  equation: 


dgjp, ^ 

dn 


4>{^dSg-  [  g{p,^~dSg  =  2T:d>{p) 
JSc+Sm  on 


(7) 


where  Sc  is  the  surface  of  the  cavity,  5m  the  interface  between  the  coating  and  the  fluid,  p 
is  a  field  point  that  is  on  the  surface  5  =  5c  +  5m,  9  is  the  source  point  that  is  also  on  5, 
=  l/|p  ~  ^  is  the  area  element  of  5  varying  the  point  q.  Once  this  equation 

is  solved  the  calculation  can  proceed  on  to  the  next  time  step  or,  with  a  companion  form  of 
this  integral  equation,  the  velocity  and  pressure  can  be  found  at  any  point  in  the  fluid. 


3  Numerical  Implementation 

In  the  numerical  model  the  surface  of  the  cavity  is  approximated  by  a  set  of  surfaces  (panels) 
each  of  which  is  obtained  by  rotating  a  straight  line  in  the  r  —  z  plane  about  the  z-axis  (see 
Figure  2).  The  cavity  is  composed  of  nj  of  these  panels.  The  bottom  boundary  of  the  fluid 
is  modeled  by  a  is  set  of  of  these  panels  located  at  z  =  0  with  normal  in  the  z-direction. 
Assume  that  at  any  time,  ti,  (the  subscript  is  used  to  denote  the  time  step)  the  dependent 
variables  describing  the  flow  and  the  membrane  wall  are  completely  know  and  it  is  desired 
to  advance  to  the  next  time  step,  L+i.  The  r,  z  coordinates  at  of  the  n;,  -f  1  intersections 
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(nodes)  of  the  lines  joining  adjacent  panels  on  the  cavity  with  the  r  —  z  plane  are  found  from 
the  finite  difference  form  of  Equations  (3): 

rUi  =  ri -  U)  (8) 

-^)  (9) 

where  the  superscript  j  refers  to  the  node  (see  Figure  2).  The  value  of  0  at  U+i  at  each  node 
is  obtained  from  the  finite  difference  form  of  Equation  (2) 

=  <»;  +  (^4.  -  <.)  (lo) 

At  this  point  we  have  the  position  of  the  nodes  on  the  cavity  and  the  value  of  0  at  each  node 
at  t,+i .  The  next  step  is  to  obtain  90  dn  at  the  nodes  on  the  membrane-flow  interface. 


The  first  step  in  the  solution  of  the  membrane  problem  is  to  Hankel  transform  Equa¬ 
tion  (4).  After  combining  Equation  (6)  wnth  the  membrane  equation  (4)  and  performing  the 
Hankel  transform,  the  membrane  equation  becomes 

Kfj  =  -P(k,0,i)  (11) 

where  for  any  quantity  Q  the  Hankel  transform  is  given  by: 

Q{k,z,i)=  /  Q{r,z,i)rJo{kr)dr  (12) 

and  k  is  called  the  wavenumber.  The  values  of  fji+i  at  each  wavenumber  k^  are  obtained 
from  the  temporal  finite  difference  form  of  Equation  (11)  : 

nL  =  (2 + + K)—)fii  -  vL,  -  (13) 

m  m 

The  values  of  at  the  nodes  between  wall  panels  are  then  obtained  by  inverting  the  Hankel 
transform.  The  Hankel  transforms  are  performed  using  a  numerical  algorithm  described  by 
Anderson  [13].  The  vertical  velocity  at  the  coating  surface  is  then 


dri  90 

dt  dn 


(14) 


Thus  after  this  computation  the  position  and  the  valujs  of  0  are  known  for  the  nodes  on 
the  cavity  and  the  values  of  90/9n  are  known  for  the  nodes  on  the  membrane.  In  order  to 
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advance  to  the  next  time  step  the  values  of  d(i>ldn  on  the  cavity  surface  and  the  values  of  <i> 
on  the  membrane  surface  must  be  computed  from  Equation  (7).  In  the  numerical  solution 
of  this  integral  equation  the  values  of  <f)  and  d<j>ldn  are  assumed  to  vary  linearly  in  the  r,  2 
plane  along  each  panel.  Thus,  in  its  discrete  form  the  integral  equation  becomes 


where  N  —  nc-¥  rim  and 


=  /  dV 

Jo 


j: 


GN 


Jj  Jo  Ip*  —  ^)| 

fV  .  li  f2r  1 

=  /  dP—  [  -z,, = - de 

Jo  Lf^  Jo  lo*  —  a3 

.0=1 


de 


dV 


Ip*  -  q\r,e)\ 

D-V  /2ir  f)  1 


LJ 


/„ 


dn  |p‘  _  q^{T^ e)\ 


dO 


...  P  d 


GN^’'’^  =  /  dl^j-  [ 
Jo  LJ  Jo 


-dO 


(16) 

(17) 

(18) 

o  -  -  (19) 

dn  |p<  _  qJ(r,0)| 

In  these  equations  the  length  of  panel  j  in  the  r  —  2  plane  is  given  by  P  and  dP  is  the 
corresponding  differential  length  element.  The  parameter  o’  is  a  geometrical  parameter.  On 
the  flow-membrane  interface,  which  is  taken  at  2  =  0,  a  tangent  plane  to  the  paneled  surface 
exists  at  each  node  and  o^  =  27r  for  these  nodes  (ric  +  2  <  j  <  -|-  n,„  +  2)  els  it  would  on  a 

continuous  surface.  However,  for  the  nodes  on  the  cavity  surface,  there  is  no  tangent  plane 
and  o'  =  lim<_oo  Aj/e^  where  e  is  the  radius  of  a  small  sphere  centered  on  the  node  i  and 
/!/  is  the  portion  of  the  surface  area  of  the  sphere  that  is  in  the  fluid. 


Starting  the  calculations  at  t  =  0  has  turned  out  to  be  rather  difficult.  The  problem 
lies  in  the  value  of  the  pressure  on  the  membrane  at  f  =  0.  Since  the  values  of  (f)  on  the 
cavity  and  d<t>ldn  on  the  membrane  are  zero  at  <  =  0,  use  of  the  integral  equation  will  yield 
^  =  0  on  the  membrane  at  this  time.  Thus,  from  Bernoulli’s  equation,  the  pressure  on  the 
membrane  at  f  =  0  is  determined  by  d(f)fdt.  To  compute  this  derivative  the  values  of  <f> 
at  f  =  At  must  be  determined.  In  proceeding  to  time  At  to  obtain  the  required  values  of 
<f>  on  the  membrane,  the  position  of  the  nodes  on  the  cavity  and  the  corresponding  values 
of  <i>  are  obtained  from  Ekjuations  and  (8)  and  (10).  In  order  to  get  4)  on  the  membrane 
surface  the  value  of  d<i>dn  must  be  known  there  so  that  the  integral  equation  (15)  can  be 
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solved.  However,  d(i>dn  is  computed  from  the  membrane  velocity  and  this  value  requires  the 
pressure  on  the  membrane  at  /  =  0.  Unfortunately,  this  pressure  is  the  quantit}-  tliat  we 
have  been  trying  to  compute  in  the  first  place.  An  iterative  scheme  was  devised  to  o\’ercome 
this  problem.  On  the  first  iteration,  the  pres'iure  on  2  =  0  at  t  =  0  is  computed  assuming  the 
coating  to  be  rigid.  This  allows  for  the  computation  of  d4>dn  on  the  membrane  at  t  =  At. 
It  hais  been  found  that  the  distribution  of  the  pressure  on  the  membrane  oscillates  from  one 
time  step  to  the  next  with  this  starting  pressure  assumption.  However,  if  the  membrane 
wall  is  stiff  enough,  the  pressure  seems  to  settle  down  to  a  constant  distribution  after  several 
time  steps.  For  these  stiff  coatings,  the  calculation  was  restarted  at  this  point  with  the 
equilibrated  pressure  distribution  rather  than  the  rigid  wall  distribution  used  at  f  =  0.  For 
walls  that  are  more  compliant,  the  pressure  distribution  diverges  when  starting  with  the 
rigid  wall  pressure  distribution.  Thus,  in  order  to  perform  the  calculation  in  these  cases, 
the  calculation  was  started  with  a  stiff  wall  and  after  several  time  steps,  when  the  pressure 
equilibrated,  the  calculation  was  restarted  with  less  stiff  wall  parameters.  This  procedure 
was  repeated  until  the  desired  wall  properties  were  reached. 


4  Results 

There  are  a  number  of  independent  variables  in  the  the  cavity  collapse  problem.  For  the 
fluid,  these  variables  include  the  initial  cavity  radius,  and  distance  from  the  wall,  Zq,  the 
pressure  difference,  Poo  —  Po  and  the  density  of  the  fluid,  p.  For  the  membrane  the  variables 
are  the  mass  per  unit  area,  m,  the  tension,  P,  and  the  spring  constant,  K.  Thus,  the  collapse 
time,  Tc,  is  a  function  of  seven  independent  parameters  which  contain  the  dimensions  of 
mass,  length  and  time.  In  dimensionless  form  there  are  five  independent  parameters  and  if 
we  define  To  =  Ro\JpI{P<x>  —  Po)i  (the  characteristic  time  scale  for  the  collapse  of  a  cavity 
in  an  infinite  fluid)  and  ~  yjmj K,  (the  characteristic  spring-mass  time  scale  for  the 
coating),  we  find: 

h- fib.  Zk  bl 

To  ^  VZo’  To’  T  ’  AP  j  ^  ^ 

where  AP  =  Poo  —  Po-  In  the  following  we  have  varied  only  KRqIAP  with  the  other  dimen¬ 
sionless  parameters  held  constant:  RqIZq  =  2/3,  Tm/To  =  1,  and  P^K fT  =  12.5.  In  these 
calculations,  the  cavity  is  composed  of  32  panels  which  are  of  equal  arc  in  the  r  —  2  plane 
when  the  cavity  is  spherical.  There  are  70  panels  on  the  flow-membrane  interface.  These 


panels  do  not  move  from  there  positions  at  <  =  0.  The  first  node  on  this  interface  is  at  r  =  0, 
and  the  next  is  at  r  =  0.05/?o-  The  length  of  subsequent  panels  increases  linearly  with  the 
laist  node  at  r  =  100.0/?o-  In  these  initial  calculations  the  time  step  was  constant  throughout 
the  collapse  and  equal  to  0.00125ro.  Because  the  rate  of  change  of  (f)  increases  rapidly  toward 
the  end  of  the  collapse  while  the  time  step  is  constant,  the  calculations  become  inaccurate 
just  before  the  collapse  is  complete.  It  is  estimated  that  inaccuracies  become  important 
at  about  a  time  of  0.02ro  before  the  end  of  the  collapse.  At  the  present  time  a  variable 
time  step  version  of  the  computer  program  is  under  development  to  remedy  this  problem. 
However,  the  results  in  the  following  discussion,  which  show  for  instance  the  collapse  time 
of  the  cavity  as  a  function  of  coating  properties,  are  still  valid  for  comparison  purposes. 

Figure  3  contains  profiles  of  the  cavity  at  various  times  during  the  collapse  for  the  ca.se  of 
a  rigid  wall,  KRo/AP  =  oo,  and  the  case  of  a  compliant  wall  KRoJAP  =  12.5.  Note  that 
the  profiles  appear  fairly  similar  except  that  the  collapse  occurs  over  a  longer  time  for  the 
case  with  the  rigid  wall.  Note  also  that  the  vertical  extent  of  the  cavity  at  the  final  instant 
is  greater  for  the  case  with  the  rigid  wall  and  that  the  height,  where  the  upper  (north) 
pole  of  the  cavity  meets  the  lower  (south)  pole  of  c.'vity  is  larger  for  the  compliant  wall  case. 
To  emphasize  this  later  point  plots  of  the  the  height  of  the  north  and  south  poles  of  the 
cavities  versus  time  for  the  rigid  wall  case  and  cases  with  KRo/AP  =  50.0, 25. 0,  and  12.5 
are  shown  in  Figure  4.  Note  the  trend  toward  larger  values  of  Zc  and  shorter  collapse  times 

>  as  KRo/AP  is  decreased  from  oo  (the  a  rigid  wall)  to  12.5.  For  the  rigid  wall,  z^  =  1.035/?o 
and  Tc  =  1.0357bi  while  for  the  case  with  KRq/AP  =  12.5  ,  z^  =  1.128i?o  and  Tc  =  0.936To. 
In  the  case  of  a  cavity  collapsing  in  an  infinite  fluid  the  collapse  time  is  0.9157^0  and  would 
equal  Zq  =  1.5.  Thus,  the  trend  is  toward  the  infinite  fluid  case  as  KRq/AP  is  decreased. 

► 

The  pressure  distributions  (P(r)  —  Poo)/{Poo  —  Po)  on  the  flow-membrane  interface  at 
times  corresponding  to  Figure  3  are  shown  in  Figure  5  for  the  case  with  the  rigid  wall  and 
)  the  case  with  KRq/AP  =  12.5.  Note  that  the  plots  only  extend  to  r  =  20i?o  while  the 

wall  pamels  extend  to  r  =  lOO/Zo-  In  both  cases  the  pressure  is  initially  slightly  negative 
directly  under  the  cavity  ajid  becomes  positive  with  a  value  of  about  15.0  toward  the  end  of 
the  collapse.  At  first  glance  the  shapes  of  these  distributions  seem  fairly  similar.  However, 
^  the  values  at  t/Rq  —  20.0  are  negative  on  the  compliant  wall  and  positive  on  the  rigid  wall. 

Since  the  horizontal  derivative  of  the  pressure  is  about  the  same,  the  component  of  the  flow 
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tangent  to  the  wall  is  outward  and  must  be  of  about  the  same  magnitude  in  the  two  cases. 
In  the  compliant  wall  case  the  gradient  must  eventually  change  sign  in  order  to  reach  zero  at 
r  =  oo  indicating  flow  toward  the  2— axis.  In  Figure  6  the  pressure  on  the  membrane  at  r  =  0 
is  plotted  versus  time  for  the  case  with  the  rigid  wall,  the  cases  with  KRq/AP  =  50.0,25.0, 
and  12.5,  and  a  case  with  '  collapse  in  an  infinite  fluid  (spherical  collapse).  In  the  latter 
case  the  pressure  wcis  evaluated  at  a  distance  Zq  from  the  center  of  the  spherical  cavity. 
Note  that  all  the  curves  have  about  the  same  shape,  but  as  the  wall  is  softened  from  a  rigid 
wall  to  KRoJAP  =  12.5  and  finally  to  the  spherical  case  (no  wall),  the  curve  is  shifted  up 
slightly  and  to  the  left.  This  is  a  further  indication  that  softening  the  wall  is  making  the 
characteristics  of  the  collapse  tend  toward  the  spherical  collapse  case. 

The  distribution  of  vertical  velocity  of  the  membrane  surface  at  various  times  is  plotted 
in  Figure  7  for  two  cases.  In  Figure  7a  the  full  membrane  coating  interaction  was  computed 
for  KRafAP  =  12.5.  However,  in  Figure  7b  the  pressure  from  the  rigid  wall  calculation 
was  used  to  drive  the  membrane,  also  with  KRqJ AP  =  12.5.  (The  membrane  motion  was 
not  ted  back  to  the  flow  model  in  this  partial  interaction  case.)  Again  the  curves  for  the 
two  cases  look  similar.  However,  the  velocities  for  the  partial  interaction  Ccise  are  larger  at 
r  =  20i2o-  The  maximum  value  of  the  membrane  displacements  corresponding  to  this  data 
is  about  O.O25R0.  Since  the  horizontal  length  scale  of  the  membrane  deformations  is  Ro,  the 
deformations  are  of  exceedingly  small  slope  and  the  linearization  of  the  boundary  conditions 
at  the  flow  coating  interface  is  justified.  Plots  of  the  velocity  at  r  =  0  versus  time  for  the  full 
interaction  cases  with  KRojAP  =  50.0,25.0,  and  12.5,  the  partial  interaction  case,  and  the 
fluid  velocity  at  a  distance  Zq  from  a  spherical  collapse  are  shown  in  Figure  8.  Note  that, 
as  expected,  the  velocity  increases  with  time  and  increases  as  the  spring  stiffness,  K RofAP, 
decreases.  The  partial  interaction  case  and  the  corresponding  full  interaction  case  are  identi¬ 
cal  until  about  t/To  =  0.6.  However,  they  are  quite  different  at  the  end  of  the  collapse.  The 
maximum  value  of  the  velocity  for  the  spherical  collapse  is  roughly  five  times  larger  than  the 
value  for  the  full  interaction  with  KRq/AF  =  12.5.  Thus,  under  the  assumption  that  the 
collapse  will  become  spherical  when  ihe  compliant  wall  behaves  in  a  manner  similar  to  the 
fluid  at  that  position  in  a  spherical  collapse,  the  wall  will  have  to  be  considerably  softer  to 
prevent  the  formation  of  the  jet  directed  toward  the  wall. 
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5  Conclusion 


A  numerical  method  for  the  computation  of  the  collapse  of  a  cavity  adjacent  to  a  compliant 
surface  has  been  presented.  The  results  of  preliminary  calculations  show  that  as  the  spring- 
backed  membrane  becomes  soft,  the  collapse  characteristics  change  in  the  direction  of  the 
collapse  of  a  spherical  cavity  in  an  infinite  fluid.  At  the  present  time  calculations  for  a 
wider  range  of  conditions  are  being  pursued  and  the  calculation  method  is  being  modified  to 
include  variable  time  steps  and  computations  of  the  energy  exchange  between  the  flow  and 
the  coating. 
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Figure  1  Schematic  showing  the  coordinate  system  and  the  initial 
position  of  the  cavity  and  the  compliant  boundary 


Figure  2  Nodes  and  panels  on  the  surface  of  the  cavity  and 
the  flow  membrane  interface. 
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Figure  4.  Height  of  North  and  South  poles  of  the  cavity  versus  time  for  four  cases,  a: 
RojZo  =  1.5,  T^ITo  =  1.0,  RlKIT  =  12.5,  KR^jAP  =  12.5.  b:  Ro/Zo  =  1.5,  T^/To  =  1.0, 
R^K/T  =  12.5,  KRo/AP  =  25.0.  c:  R^fZ^  =  1.5,  TJTo  =  1.0,  RlK /T  =  12.5, 
KRofAP  =  50.0.  d:  rigid  wall,  Rq/Zq  =  1.5. 


17 


Figure  6.  Pressure  at  r  =  0  on  the  membrane  surface  versus  time,  a:  the  pressure  at  the 
radius  Zq  for  a  spherical  cavity  collapsing  in  an  infinite  fluid,  b;  /?o/^o  =  1-5,  T^ITo  =  1.0, 
n^K/T  =  12.5,  KI^I^P  =  12.5.  c;  B^jZo  =  1.5,  r,„/To  =  1.0,  RlKlT  =  12.5, 
KR^IAP  =  25.0.  d:  RofZo  =  1.5,  T^/To  =  1.0,  P^KfT  =  12.5,  KR^IAP  =  50.0.  e: 
rigid  wall,  RoIZq  =  1.5. 
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locity  at  the  radius  Zq  for  a  spherical  cavity  collapsing  in  an  infinite  fluid,  b:  Rq/Zq  =  1 
TJTo  =  1.0,  R^K/T  =  12.5,  KR^^IAP  =  12..5.  c:  R^fZo  =  1.5,  Tm/To  =  1.0,  R^K/T 
12.5,  KRo/AP  =  25.0.  d;  R^/Zq  =  1.-5,  T„,/To  =  1.0,  RlK/T  =  12.5,  KR^/AP  =  50.0. 
rigid  wall,  R^lZo  =  1.5. 


